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Abstract 

We investigate phenomenology related to the neutral fields in supersymmetric models 
with an extra U(l) derived from E e . Our study is concentrated into the models which 
have a singlino dominated neutralino as the lightest superparticle (LSP). If such models 
satisfy a constraint for dark matter derived from the WMAP data, the lightest neutral 
Higgs scalar, a new neutral gauge field Z' and the LSP may be interesting targets for the 
study at the LHC. We also discuss features of the Z' in the models and its detectability 
at the LHC. 
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1 Introduction 



Recent various astrophysical observations quantitatively show the existence of a substan- 
tial amount of non-relativistic and non-baryonic dark matter [1, 2]. Although super- 
symmetric models have been considered to be the best candidate beyond the standard 
model (SM) from a viewpoint of both the gauge hierarchy problem and the gauge coupling 
unification, this fact seems to make them much more promising on the basis of an ex- 
perimental signature [3]. If -R-parity is conserved in supersymmetric models, the lightest 
superparticle (LSP) is stable. Thus the LSP can be a good candidate for cold dark matter 
(CDM) as long as it is electrically neutral. Since the strength of interactions of the LSP 
with the SM fields is 0(Gp) and the mass can be of the order of the weak scale, its relic 
energy density is expected to be eventually of the order of critical energy density of the 
universe. The most promising one among such LSP candidates is considered to be the 
lightest neutralino. Relic abundance of the lightest neutralino has been extensively stud- 
ied in the minimal supersymmetric SM (MSSM) [4, 5, 6]. After publication of the analysis 
of the WMAP data, however, the allowed parameter space in the minimal SUGRA (or 
CMSSM) is found to be restricted into some narrow regions [7]. If we take account of this 
situation, it seems worth studying the relic abundance of CDM candidates quantitatively 
also in various extensions of the MSSM. It may also be useful to discuss indications of 
such models which are expected to be found at the forthcoming Large Hadron Collider 
(LHC) by applying the CDM condition. 

The MSSM has been considered as the most promising supersymmetric model and 
has extensively studied from various points of view. Although the MSSM can explain 
experimental results obtained by now as long as parameters are suitably chosen, it suffers 
from the well known \i problem [8]. If we try to solve it near the weak scale, we need to 
extend the MSSM in the way to give some influence to physics at TeV regions [9, 10]. If 
we extend it by introducing an extra U(l) gauge symmetry, for example, the problem can 
be solved in a very elegant way [10, 11]. The existence of additional U(l) symmetries is 
also predicted in some effective theories of superstring [12, 13]. If there is an extra U(l) 
symmetry at TeV regions, the model is expected to reveal distinguishable features from 
those of the MSSM. Some signals of the models may be detected at the LHC [13, 14, 15] 
and a CDM candidate may be different from that in other models. In this paper we focus 
our attention into such aspects in the models with an extra U(l). 
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In the models with an extra U(l), an operator \SHiH 2 is introduced into superpo- 
tential as a gauge invariant operator instead of the so-called pi term [iHiH 2 [10]. 1 In the 
simplest models to accommodate such a feature, the extra U(l) symmetry is supposed 
to be broken by a vacuum expectation value (VEV) of the scalar component S of an SM 
singlet chiral superfield S. The \i term is generated as ji = X(S) by the same singlet scalar 
field through the introduced operator. These models show a difference from the MSSM in 
the neural Higgs sector in addition to the existence of a new neutral gauge field Z' [10, 16]. 
Since the neutralino sector is extended from the MSSM by a fermionic component S and 
an extra U(l) gaugino X x , the feature of neutralinos can also be different from that in the 
MSSM. Various phenomena are influenced by this change [17]. 

In particular, if the singlino S can dominate the lightest neutralino, distinguishable 
features from the MSSM are expected to appear in the phenomena relevant to the neu- 
tralinos. When (S) takes a value of the weak scale, the singlino domination of the lightest 
neutralino is naturally expected to occur. This is the case in the well known next MSSM 
(NMSSM) [18] and its modified model (nMSSM) [19]. In the models with an extra U(l), 
however, there are severe constraints on (S) from both mass bounds of Z 1 which result 
from the direct search of Z 1 [20] and bounds on the mixing between Z 1 and the ordinary Z 
which result from the electroweak precision measurements [13, 21, 22]. These constraints 
tend to require that (S) should be more than 0(1) TeV as long as we do not consider 
a special situation. 2 Thus, in the simple models with an extra U(l) which is called the 
UMSSM in [24], it seems unable to expect substantial differences in the lightest neutralino 
from the MSSM since both S and tend to decouple from the lightest neutralino. 

This situation changes if the extra U(l) gaugino X x is sufficiently heavier than Z' as 
suggested in [11, 23, 24, 25]. In these papers the nature of the lightest neutralino has 
been studied to show a lot of interesting aspects: (a) it can be dominated by the singlino 
S and can be very light; (b) it can be a good candidate of the CDM whose nature is 

1 In this paper we put a hat on the character for a superfield. For its component fields, we put a tilde 
on the same character to represent the superpartners of the SM fields and use just the same character 
without the hat for the SM fields. Otherwise, the field with no tilde should be understood as a scalar 
component. 

2 Even in the models with an extra U(l), if one considers a model with a secluded singlet sector which 
is called the S'-model in [24], (S) can take a value of the weak scale. In this case phenomcnological feature 
at the weak scale is very similar to the nMSSM. 
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very different from that in the MSSM. It can have a small mass which has already been 
forbidden in the MSSM since it is dominated by the singlino; (c) it also has a different 
interaction from that in the MSSM. However, the studied parameter regions are different 
in each study. In [23] the CDM abundance is mainly studied under the assumption that 
the lightest neutral Higgs mass is m h = 170 GeV. In the S-model case the neutralino mass 
matrix is assumed to be reduced to the one of the nMSSM [24] or a certain gaugino mass 
relation such as M x = M Y is assumed [26] . In this paper we do not use these assumptions. 
We are interested in other region of the parameter space, where X x is much heavier than 
other gauginos but it has non-negligible mixing with S [25, 27]. Therefore \ x can give a 
crucial contribution to the mass of S. We focus our study on such cases and investigate 
the neutral field sector in the models with an extra U(l) derived from E % . We discuss 
phenomenological features of the models and study indications of the models expected to 
be found at the LHC. 

The paper is organized as follows. In section 2 we define the models and discuss their 
features different from other models. In particular, we focus our attention on the neutral 
field sectors of the models, that is, the lightest neutralino, the new neutral gauge field Z' 
and the neutral Higgs scalars. In section 3 we study the parameter space of the models 
which is allowed by the current results of various experiments. We predict signals of the 
models expected to be seen at the LHC. Section 4 is devoted to the summary. 

2 11 problem solvable extra U(l) models 
2.1 Features of the neutral fields 

In the simple models with an extra U(l) which can solve the /i problem, the mass of Z' is 
directly related to the /x term [10, 11]. This feature induces various interesting phenomena 
which make the models distinguishable from the MSSM. The [i term is considered to be 
generated by an operator in the last term of the superpotential 

W oh = huUQH 2 + hoDQHi + h E ELH x + \SH X H 2 , (1) 

where are the ordinary doublet Higgs chiral superfields and S is an additional singlet 
chiral superfield. If a VEV of the scalar component of S is assumed to generate both 
the /i term and the Z' mass, the superpotential (1) requires that Hi 2 also have the extra 
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U(l) charges Qi t2 . The charge conservation imposes a condition Q\ + Q2 + Qs = on 
them. A bare \i term nHiH 2 is automatically forbidden by this symmetry. Stability of 
the scalar potential for S is automatically guaranteed by a quartic term induced as the 
extra U(l) D-term without introducing a new term in the superpotential. The Z' mass 
is difficult to be much larger than 0(1) TeV as long as there is no other contribution 
to it. These aspects may make the models not only theoretically interesting but also a 
good target for the studies at the LHC [13, 14, 15]. In particular, we can find interesting 
features in the neutral fields, which make the models distinguishable from the MSSM 
and also the singlet extensions of the MSSM such as the NMSSM, the nMSSM and the 
S-model. In this section we review these features to clarify the differences of our models 
from previously studied ones. 

Before proceeding to the discussion on this issue, we fix assumptions for the super- 
symmetry breaking. We assume soft supersymmetry breaking terms 

-£susy = E m jM 2 + [\ M Ag~K + \m w \ w \ w + ^My\ Y \ Y + ^M x \ x \ x + h.c.) 

(AuhuUQH 2 + A D h D DQHi + A E h E ELH x + A X \SH X H 2 + h.c.) , (2) 

where tp in the first term runs all scalar fields contained in the models. Scalar mass m 9 
and A parameters for scalar trilinear terms may be assumed to have a universal value 
TO3/2. A grand unification relation for the masses of gauginos may also be assumed. 

Now we assume that the scalar components of Hi j2 and S obtain the VEVs t>i, v 2 and 
u due to radiative effects on the supersymmetry breaking parameters and the electroweak 
symmetry SU(2) L xU(l)y xU(l) x breaks down into U(l) cm [10, 11]. We define the field 
fluctuations around this vacuum as 

" L ' ,B1 ' " 1 I , H 2 =\ 7 V S = u + h° s + iP s . (3) 





The last term in eq. (1) generates a \x parameter as fi = Xu and plays a required role for the 
\i term in the MSSM. Although the models have similar structures to the MSSM after this 
symmetry breaking, there appear various differences induced by this term, in particular, 
in the sector of the neutral fields. In a charged field sector, we can find a difference from 
the MSSM only in the mass of the charged Higgs scalar. It will be discussed in Appendix. 
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2.1.1 Neutral gauge field sector 

The clearest difference from the MSSM is the existence of the Z' at TeV regions. Through 
the symmetry breaking denoted by eq. (3), a mass matrix M\ z , is generated for the neutral 
gauge bosons Z^ and Z^. It can be written in the basis (Z^, Z^) as 
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-v 2 (Q lC os 2 (3-Q 2 sm 2 (3) 



(4) 



V 9xV Y ?1 v 2 (Q 1 cos 2 p-Q 2 sm 2 p) §v 2 (Q 2 cos 2 /? + Q 2 sin 2 (3 + Q| J) 

where i> 2 = t> \ + tf and tan/3 = v 2 /v±. The extra U(l) charge Q x (f) of a field / and its 
coupling g x are defined in such a way that the covariant derivative for the relevant gauge 
group SU(2) L xU(l)y xU(l) x takes the form 

D, = d,- igJ-^W; - ig^f-B, - ig.^ 1 ^. (5) 

Mass eigenvalues of the neutral gauge fields Z lfJi and Z 2li can be approximately expressed 
as 

~ m z — m z —j^^^={Qi cos (3 - Q 2 sm p), 

\f92 + gi 

< ^ 9 4(QWi + Qlvl + QW) + m|^%Qi cos 2 /? - Q 2 sin 2 /?), (6) 

2 \Jgl + 9i 

where mz is the mass of the Z boson in the SM and £ is a ZZ' mixing angle defined by 

= 2^ V / ( ?2 2 + ( ?1 2 (Q 1 cos 2 /3-Q 2 sin 2 /j) 

^ 4 <£(Q? cos 2 (3 + Q 2 2 sin 2 /3 + Q 2 s u 2 /v 2 ) - (g 2 + g 2 ) ' {V 

If -u ^> v i 5 2 is satisfied, approaches to and mz 2 is proportional to u. The relation 
between m^ 2 and u is given by mz 2 /u ~ g x Qs/V^- 

Both direct search of a new neutral gauge field and precise measurements of the elec- 
troweak interaction severely constrain the mass eigenvalue mz 2 of the new neutral gauge 
boson and the ZZ' mixing angle £ [21]. Lower bounds for m Z2 have been studied by using 
the searches of the Z 2 decay into dilepton pairs [20] . Although it depends on the models, 
it may be roughly estimated asm Zj > 600 GeV. If Z 2 has a substantial decay width 
into non-SM fermion pairs such as neutralino pairs, this bound may be largely relaxed 
[15]. On the other hand, the precise measurements of the electroweak interaction give a 
constraint £ < 10~ 3 [21]. As found from eq. (7), this bound can be fulfilled if either of 
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two conditions is satisfied, that is, a sufficiently large u or tan/? ~ \JQ\IQ2 [H]- For the 
latter case, since the constraint from the ZZ' mixing can automatically be guaranteed, u 
needs not so large as long as the direct search constraint on m Z2 is satisfied. This seems 
to be an important point to be noted when we consider the existence of an extra U(l) at 
TeV regions. We will focus our study on this case, which may make it possible to find 
the solutions in the parameter regions excluded in the study of the UMSSM [23, 24]. 

2.1.2 Neutralino sector 

The neutralino sector is extended into six components, since there are two additional 
neutral fermions X x and S. X x is the extra U(l) gaugino and S is the fermionic component 
of S. If we take a basis Af T = {—iX x , —iX^, —iX Y) H 1: H 2 , S) and define neutralino mass 
terms such as £^ utralino = — \N T M.N + h.c, a 6 x 6 neutralino mass matrix M. can be 
represented as 3 

( M~ x ^vcos/3 ^vsin/5 9 -^u \ 

Mjy Tn z cw cos (3 —m z cw sin (3 

My —rnzsw cos (3 m z sw sin (5 

a ^-v cos P rn z cw cos (3 —m z sw cos (3 \u Xv sin (3 

^%-vsiiif3 —m z cw sin (3 m z sw^(3 Xu Xvcos(3 

y a^u Xv sin (3 Xvcosfi j 

(8) 

Neutralino mass eigenstates Xa(a = 1 ~ 6) are related to A/} by using the mixing matrix 

U as 

X° a = T,U aj Af„ (9) 
j'=i 

where U is defined in such a way that UMU T becomes diagonal. 

The composition of the lightest neutralino is important for the study of various phe- 
nomena, in particular, the relic density of the lightest neutralino as a CDM candidate. 
If u is a similar order value to t>i 2 or less than these, the lightest neutralino is expected 
to be dominated by the singlino S just like in the case of the nMSSM and the S-model. 
In this case, if it can annihilate sufficiently, the lightest neutralino with a sizable singlino 



3 We do not consider gauge kinetic term mixing between U(l)y and the extra U(l), for simplicity. The 
study of their phenomenological effects can be found in [17]. 
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component may be a good CDM candidate in the parameter regions different from the 
ones in the MSSM [18, 19]. In the present models, however, the Z' constraints seem to 
require that u should be much larger than as mentioned before. As the result, A^ and 
S tend to decouple from the lightest neutralino as long as the mass of \ x is assumed to 
be similar to the masses of other gauginos. In such a situation, the composition of the 
lightest neutralino is expected to be similar to that of the MSSM. Then we cannot find 
distinctive features in the lightest neutralino. However, as suggested in [11, 23, 24, 25], 
this situation can be drastically changed if the mass of the extra U(l) gaugino M x be- 
comes much larger than the masses of other gauginos due to some reasons. In this case 
the lightest neutralino can be dominated by the singlino S. 

If the gaugino X x is heavy enough to satisfy M x 3> ^^m, we can integrate out A^ 
just as the seesaw mechanism for the neutrinos. A resulting 5x5 mass matrix can be 
expressed as 











m z C]Y cos (3 —m z sw cos (3 
—m z cw sin (3 m z sw sin (3 



\ 











m z cw cos (3 
—m z sw cos (3 

Am 
Xv sin (3 



—m z cw sin (3 
m z sw sin (3 

Am Xv sin (3 

-f^Jv 2 sin 2 Xvcos(3 

- 9 M u2 



\ 



(10) 



This effective mass matrix suggests that the lightest neutralino tends to be dominated 

~ 2q2 u 2 

by the singlino S as long as M^ y an d M= i s n °t smaller than ^m- • Since 
and u cannot to be less than 100 GeV because of mass bounds of the lightest chargino 
and the gluino [28], the singlino domination of the lightest neutralino is expected to be 
naturally realized in the case that M x ^> u is satisfied. In such a case, phenomenology of 
the lightest neutralino can be largely changed from that in the MSSM and also its singlet 
extensions. This effective mass matrix reduces to the one of the nMSSM in the large M x 
limit [24]. However, we are interested in the intermediate situation where the effectively 
generated diagonal elements in eq. (10) can not be neglected. Since this possibility has 
not been studied in detail in realistic models yet despite it is potentially interesting, 4 we 
will concentrate our study into such a situation in this paper. The lightest neutralino is 
assumed to be dominated by the singlino because of a large M x compared with the mass 



4 In particular, the detailed study seems to have not been done under the assumption tan (3 = \J Q\jQi- 
In [25] we do not consider the anomaly problem of U(l) x seriously and only study a toy model. 
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of other gauginos. 



2.1.3 Neutral Higgs scalar sector 

The neutral Higgs mass is also modified from that in the MSSM as in the case of the 
NMSSM and the nMSSM. A difference from the latter ones is the existence of an additional 
D-term contribution of the extra U(l). The CP even neutral Higgs sector is composed of 
the three scalars (h^.h^h^) which are introduced in eq. (3). Their mass matrix M.\ at 
tree level is written as 

( \{gl+gl+g 2 x Ql)vl+A x \ut&n(i -\{g\ + g 2 - g 2 (l, 2)) Vl v 2 - A x Xu ± 9 2 (1, S)v lU - A x Xv 2 \ 

-\{gi+gl-g 2 {l^)) Vl v 2 ~ A x \u \{gl + g\ + g 2 x Ql)vl + A x \ucotp \g 2 {2, S)v 2 u - A x \ Vl , 
V y 2 (l,S)v lU -A x \v 2 \g 2 (2,S)v 2 u-A x \v 1 \g 2 x Q% u 2 + A X X^ J 

(11) 

where A x is the soft supersymmetry breaking parameter defined in eq. (2). We use a 
definition g 2 (i,j) = glQiQj + 4A 2 in this formula. Mass eigenstates <p a are related to the 
original neutral CP even Higgs scalars h Q a by 

a = E Oaahl, (12) 
a=l,2,S 

where the orthogonal matrix O is defined so as to diagonalize the neutral Higgs mass 



02 ' m <p3 



matrix M h in such a way as OM h O = diag(ml,m 

Since upper bounds of the mass eigenvalue for the lightest neutral Higgs scalar h° can 
be estimated by using eq.(ll) as 



2 ^ 2 

m h o < m z 



2A 2 a 2 

cos 2 2(3 + -„ , sin 2 2(3 + J* x A Q X cos 2 (3 + Q 2 sin 2 (3f 

92+9i 92 + 9i 



+ Am 2 , (13) 



it can be larger than that in the MSSM. The second term in the brackets of eq. (13) comes 
from the interaction given by the last term in eq. (1). It can give a large contribution for 
smaller values of tan (3 and u for a fixed fi. The third term is the D-term contribution 
of the extra U(l). Due to these effects, even in the regions of the small tan (3 such as 
tan/3 = 1 — 2, the mass of the lightest neutral Higgs m$ a can take larger values than 
that in the MSSM, such as 120 GeV or more, if one-loop corrections Am 2 are taken into 
account [16]. Since dominant components of this lightest Higgs scalar are expected to 
be h\ 2 as long as u is not smaller than f lj2 , its nature is similar to that in the MSSM, 
except that it is heavier than that in the MSSM. This situation seems very different from 
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the nMSSM and the S-model where the lightest neutral Higgs is a mixture state of h° 12 
and h° s . In numerical studies given in section 3, we will estimate both the Higgs mass 
eigenvalues and their eigenstates by diagonalizing the mass matrix (11) including the 
one-loop corrections due to the stops. 

CP odd Higgs scalars are also somewhat changed from the ones in the MSSM. A CP 
odd Higgs mass matrix M. 2 P can be written as 



M 2 P 



( 



V 



A\Awtan/3 A\\u A\\v2 
A\\u A\\u cot 3 A\\vi 

A\\V2 A X \V! 



\ 



A X X^ 



(14) 



Only one component Pa has a non-zero mass eigenvalue 



m P A = 



2A\\u v 2 . 2 N 

7T 1 H n sin 23 

sin 2(3 1 Au 2 H ' 



(15) 



and others are would-be Goldstone bosons G® 2 as in the MSSM. This requires \uA\ > 
for the stability of the vacuum. Imaginary parts of the original Higgs fields in eq. (3) have 
Pa as a component. They can be written as 



u sin 3 
Pi = -^Pa + ■■■, 



u cos /3 „ v sin 3 cos 3 „ , . 

H -P A + . . . , P s = ^ -Pa + .-., (16) 



N N N 

where a normalization factor iV is defined as iV = \ju 2 + v 2 sin 2 3 cos 2 3. Although m 2 p A 
takes larger values than those in the MSSM, Pa is found to be similar to that in the 
MSSM if u becomes larger than v. 



2.2 A CDM constraint 

In the models with a large M £ the composition and the interaction of the lightest neu- 
tralino can be very different from that in the MSSM and its singlet extensions. Since the 
lightest neutralino can also be a CDM candidate in such models, the relic abundance is 
expected to give a different constraint on the parameter space from that in the MSSM 
and its singlet extensions. It is useful to discuss this point briefly here. 

The relic abundance of the stable lightest neutralino x° which is thermally produced 
can be evaluated as thermal abundance at its freeze-out temperature Tp. This tempera- 
ture can be determined by if (TV) ~ T^o. H(T F ) is the Hubble parameter at TV [29]. T^o 
is an annihilation rate of Xe an d it can be expressed as T^> = (a mm v )n^o, where {a ami v ) is 
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thermal average of the product of an annihilation cross section er ann and relative velocity 
v of annihilating x°s. Thermal number density of non-relativistic Xe & t this temperature 
is expressed by n^o. If we use the non-relativistic expansion for the annihilation cross 
section such as cr ami v ~ a + bv 2 and introduce a dimensionless parameter xp = m^o/Tp, 
we find that xp can be represented as 

0.0955m p im^o(a + 66/ xp) 
x F = In l —— , (17) 

{g*XFy/ 2 

where m p i is the Planck mass and enumerates the degrees of freedom of relativistic 
particles at Tp. Using this xf, the present abundance of x° can be estimated as 



oa 2 | = Xe Xe 



8.76 x 10~ u g* 1/2 x F 



(18) 



(a + 3b/ x F ) GeV 

We can find formulas of the coefficients a and b for the processes mediated by the exchange 
of various fields contained in the MSSM in the articles [4, 5]. 

If the lightest neutralino Xe is dominated by the singlino and also relatively light, the 
decay modes into other final states than the SM fermion-antifermion pairs are expected 
to be suppressed. The singlino dominated neutralino contains the MSSM higgsinos (A/4,5) 
and gauginos (A/2,3) as its components with an extremely small ratio. Thus the annihi- 
lation process caused by these components are heavily suppressed and then cannot be 
dominant modes unless the enhancement due to the pole effects of the intermidiate fields. 
As long as we consider this lightest neutralino is relatively light, such pole enhance- 
ments are kinematically forbidden in the annihilation modes which have gauge bosons 
(W + W~, Z1Z1) and Higgs scalars as the final states. Then we can expect that the anni- 
hilation through the modes XeXe ~^ ff * s dominant, since the above mentioned situation 
is escapable for this process. Moreover, exotic fields are considered to be heavy enough 
and then cannot be the final states kinematically. Thus / is expected to be restricted to 
quarks and leptons. In the present analyses we will mainly focus our attention to this 
case. These annihilation processes of the lightest neutralino in the models with an extra 
U(l) are expected to be mediated by the exchange of Zi, Z 2 and the neutral Higgs scalars 
in the s-channel and by the sfermion exchange in the t-channel. New interactions related 
to these annihilation processes of the lightest neutralino x°i can De written as 

9xQj J - / 

j=A 



E ^A>,-"A}z;, + A + « + /13A/4A/5) , 
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+ 



9xQ(f) 



+ 



9xQj 



+ 



9xQ(f) 
V2 



UlXjWvflZS + XU m (U u O a2 + U e5 O al )x° e Xe<f>a 



r-,0 7ii 



UaX°e 



//*) + •••, 



(19) 



where we use eqs. (9) and (12). 

If we consider the case that the singlino dominates the lightest neutralino, the annihila- 
tion cross section into the final states // is expected to obtain the dominant contributions 
from the exchange of the new neutral gauge field Z 2 and the exchange of the lightest neu- 
tral Higgs scalar <p a . They crucially depend on both the composition and the mass of 
and <f> a . These contributions to a and b can be expressed as [23] 



a f = 



+ 7T~ 



+ 7T~ 




1 - 



1/2 



(20) 



where Cf — 1 for leptons and 3 for quarks. The extra U(l) charges of the fermions 
are denoted by Q(fi) and Q^fn)- 'Pf ls defined by using eq. (16) as 

1 -U m O a2 (U u O a2 + U e5 O al ) (f with T 3 = i) 

(21) 



sin (3 

cos j3 



U m O al (U u O al + U e5 O a2 ) (/withT 3 = -i) 



where T 3 is the weak isospin. Contributions due to the CP odd and heavier CP even Higgs 
scalars are represented by the ellipses in eq. (20). They are expected to be suppressed 
because of their large masses. 

Since the second term of bf has no suppression from the masses of the final state 
fermions, all quarks and leptons can contribute to this term as long as the threshold is 
opened. This contribution can be effective for a larger m^o(< mz 2 ) even in the case that 
af is suppressed by a large value of m Z2 . The third term of bf is suppressed by the final 
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state fermion mass but it can give a large contribution in the case of 4> a ~ 2m^o, as 
is well known. In the present models the lightest neutralino can be much lighter than 
that of the MSSM. On the other hand, the lightest CP even neutral Higgs scalar can be 
heavier than that of the MSSM as discussed in the previous part. These features may 
make the Higgs-pole enhancement effective in the annihilation of the lightest neutralino. 
This is not realized in the nMSSM where the lightest neutral Higgs can be very light 
due to the singlet-doublet mixture. Since the lightest neutralino can have substantial 
components which have a new interaction with the lightest neutral Higgs scalar as shown 
in eq. (19), the Higgs exchange is expected to be important if the lightest neutralino 
contains the ordinary Higgsino or bino component sufficiently These aspects have been 
shown partially by using numerical studies in [25]. In addition to these effects we also have 
to take account of all other processes mediated by the exchange of the MSSM contents 
in the numerical estimation of the relic abundance of the lightest neutralino because of 
the following reasons. Firstly, the lightest neutralino with sufficient Higgsino components 
may be light enough. In that case, if m^o ~ m Zl /2 and | 14 1 2 ^ \Ui 5 \ 2 are satisfied 
as it happens in the S-model [24], the annihilation can be enhanced. Secondly, if the 
D-term contribution of the extra U(l) makes the masses of sfermions small enough, the 
t-channel exchange of those sfermions can be a crucial process for the annihilation of the 
lightest neutralinos. The D-term contribution to the sfermion mass is given in eq. (34) of 
Appendix. 

2.3 Z' decay 

Search of Z' is one of important subjects planed at the LHC except for the search of 
Higgs scalars and superpartners [13, 14]. Since the present models are characterized by 
the existence of both the new neutral gauge boson Z' and the neutral Higgs scalar heavier 
than that in the MSSM, 5 the combined analyses of these may give a useful clue for the 
search of this type of models. Here we present some useful formulas for the study of the 
Z' at the LHC. 

5 It is interesting that in a completely different context there exist other models which predict both a 
new neutral gauge boson and a neutral Higgs scalar heavier than that expected in the MSSM [30]. 
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A tree level cross section for the process pp(pp) — > Z 2 X — > ffX is given as [14] 

rl rl 

a f = J2 J q dx i J o dx 2 a(sx 1 x 2 ;qq -> ff)G q A (x 1 ,x 2 ,m 2 Z2 )0(x 1 x 2 s - M s ), (22) 

where 2^2 is defined by x = \[^-e- v using the rapidity y and a squared momentum transfer 
Q 2 . The sum of the masses of final state particles is represented by M^. s is a square 
of the center of mass energy in a collision. A function G q A (xi,x 2 ,m 2 Z2 ) depends on the 
structure functions of quarks. In the present case, eq. (22) can be approximated as [14] 



at = * 4 4^B(ff) 



B(uu) + -^-B(dd) 

G u d 



Cexp^-A^-j, (23) 



s 3 m Z2 

where C ud = 2(25), C = 600(300) and A = 32(20) for pp(pp) collisions. The QCD 
correction is taken into account by k and it is fixed to be k ~ 1.3 in the following 
numerical calculations. Yz 2 is a total width of Z 2 and B(ff) is a branching ratio of 
the Z 2 decay into //. Formulas for possible decay modes of the Z 2 are summarized in 
Appendix. 

We should note that may be expected to take rather different values from the or- 
dinary ones if the singlino dominated lightest neutralino can explain the observed CDM 
abundance. Since the decay width into the neutralino sector can be enhanced in compar- 
ison with the ordinary Z' models, the detectability of the Z 2 at the LHC may receive a 
large influence as long as the Z 2 is searched by using dilepton events (/ = e,/x) [23, 24]. 
We will compare it with the results obtained in the ordinary Z' models by practicing 
numerical analyses in the next section. 



3 Numerical analyses 
3.1 Set up for the analyses 

In this section we study parameter space of the models allowed by various phenomenolog- 
ical constraints including the CDM condition obtained from the analysis of the WMAP 
data. Then for such parameters we give some predictions of the models for the masses of 
the new neutral gauge field and the lightest neutral Higgs scalar and the detectability of 
the Z 2 at the LHC and so on. 

Before proceeding to the results of the analyses, we summarize the assumptions which 
we make in numerical studies. Firstly, we focus our study on the case that the extra U(l) 

14 
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Tabic 1 Abelian charges of the relevant fields in a fundamental representation 27 of Eq. 

gaugino has a larger mass M £ compared with other gauginos. This tends to make 
the lightest neutralino dominated by the singlino. Secondly, we consider a special case 
such that tan/? ~ \jQi/Q2 is satisfied so as to relax the constraint on the value of u 
and make a phenomenological role of the Z 2 more effective. In this case the ZZ' mixing 
constraint disappears and only the constraint derived from the direct search of the Z 2 
should be taken into account. Thirdly, we restrict our study to the extra U(l)s derived 
from Eq. The models with an extra U(l) are constrained from anomaly free conditions to 
be realistic. They are generally required to introduce the exotic matter fields to cancel the 
anomaly. If we adopt Eq as a background symmetry, we can control the field contents and 
their 11(1)^ charge systematically to make the model anomaly free. It is also a promising 
candidate for the extra U(l) models since superstring may realize them as its low energy 
effective theory [12]. 

As is well known, Eq has two Abelian factor groups in addition to the usual SM gauge 
group. We assume that only one of them remains unbroken at TeV regions and it is 
broken by the VEV u of the SM singlet scalar given in eq. (3). In this case the general 
U(l) x can be expressed as a linear combination of two representative U(l)s such as [14] 

Qx = Qip cos 6 - Q x sin 9, (24) 

where and Q x are the charges of U(l)^ and U(l) x which are obtained as 

E 6 D 50(10) x U(l)^, SO(10) d SU(5) x U(l) x . (25) 

Each charge of the relevant MSSM fields contained in the fundamental representation 27 
of E 6 is given in Table 1. The charge of other chiral superfields in the MSSM can be 
determined from them by requiring that the superpotential (1) should be invariant under 



15 




-1.5 -1 -0.5 0.5 1 1.5 



e 

Fig. 1 tan (3 for various values of 9. 

these. Although there are exotic fields such as 3 + 3* of SU(3) and 2 + 2* of SU(2) in 27 
which are not included in the superpotential (1), they can be assumed to be sufficiently 
heavy due to some large VEVs or the soft supersymmetry breaking effects. Thus, we 
neglect these effects on the annihilation of the lightest neutralino and the Z 2 decay in the 
present analyses. Since we consider that U(l) a is derived from E 6 , the coupling constant 
of U(l) x may be related to that of the weak hypercharge by g x = y§<7i, which is derived 
from the unification relation 

\gl e yf = A e e 

° /e27 /e27 
We adopt this relation in the present studies. 

In Fig. 1 we plot tan/3 for the angle 9 which is used in eq. (24) to define U(l) x . Since 
we consider the case of tan/3 = \fQ1jQ2, Q1/Q2 > 1 should be satisfied and then the 
angle 9 is found to be confined into the regions such as —0.9 ^ 9 < 0. It is interesting 
that 9 = — tan _1 (l/\/l5) — —0.253 is included in this region. Since the right-handed 
neutrinos do not have the U(l) x charge in this case, they can be very heavy without 
breaking U(l) a and the seesaw mechanism can work to realize the small neutrino mass 
[10, 22]. We will study this case in detail as an interesting example in the following. 

Now we list up free parameters in the models for the numerical analyses. In relation 
to the extra U(l) x we have the VEV u and the angle 9. Since tan f3 is assumed to be fixed 
by the U(l)a; charge as mentioned above, only A(= n/u) is a free parameter in the MSSM 
sector. For the soft supersymmetry breaking parameters, we assume the universality such 
as = m and Au — Ad = • • • = A for the parameters in eq. (2). It can be considered as 
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a result of E e . Thus, if we assume the unification relation for the masses of gauginos Mg, 
Myy, My which can be written as M g = g\M^jg\, M Y = bg\M^/Zg1, and also impose 
m = A = m 3 / 2 only to simplify the analyses, the number of remaining parameters is 
reduced into six: 6 

9, u, A, M w , M £ , 7n 3 / 2 . 

We practice the following numerical analyses by scanning the parameters u, M £ and A in 
the following regions: 

300 GeV < u < 2300 GeV (2 GeV), 200 GeV < M- x < 12 TeV (20 GeV), 
200 GeV < M w , /i < 1300 GeV (20 GeV), (26) 

where it should be noted that \i stands for A for a fixed u. In the parentheses we show 
search intervals for these parameters. We fix the supersymmetry breaking parameter m 3 / 2 
to be 1 TeV as its typical value. 7 

Throughout the analysis we impose the constraints on the masses of the chargino , 
sfermions /, the lightest CP even neutral Higgs scalar h° and the charged Higgs scalar 
/i ± as follows: 

m^± > 104 GeV, M g > 195 GeV, m f ~ > 250 GeV, 

m h o > 114 GeV, m h ± > 79 GeV. (27) 

Although the constraint on the sfermion masses are model dependent, we use the bound 
for the constrained MSSM here as an example. The sfermion masses should be checked 
whether the above bounds are satisfied by including a D-term contribution of the extra 
U(l). Since the D-term contribution may take a large negative value, this condition could 
give upper bounds for the value of u in that case. The bounds for mz 2 derived from the 
direct search of the Z 2 are taken into account by imposing the Tevatron constraint 

a(pp -> Z 2 X)B(Z 2 -> e + e-, ^ + ^) < 0.04 pb (28) 

6 It should be noted that this reduction of the model parameters are caused by the additional assump- 
tions such as tan = WQ1/Q2 and mo = A which are not related to Eq . 

7 Since is assumed to take large values, we need to consider its renormalization group effects on 

other soft supersymmetry breaking parameters and also the relation to the radiative symmetry breaking. 
However, since it depends on supersymmetry breaking scenarios, we do not go further this subject here. 
A relevant study can be found in [27] . 
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at y/s = 1.8 TeV [20]. We also impose < A = pi/u < 0.75 which is required by 
perturbative bounds for the coupling constant A [19, 24, 16]. A is also restricted indirectly 
through the relation to /j, by the chargino mass bound. The CDM constraint from the 1st 
year WMAP results is taken into account as [1] 

0.0945 < ttcDMh 2 < 0.1287 (at 2a). (29) 

If there exist other dark matter candidates than the lightest neutralino, the lower bound 
of (29) cannot be imposed to constrain the parameter region. In this paper, however, we 
make the analysis under the assumption that the lightest neutralino is the only candidate 
for the dark matter. 

3.2 Predictions of the models 

At first, in order to see global features of the extra U(l) models derived from E 6 we study 
the 9 dependence of important quantities. We scan the parameters in the regions shown 
in (26) for each value of 9 which is allowed from Fig. 1 and search the solutions which 
satisfy all conditions (27)-(29). In Fig. 2 we plot important physical quantities obtained 
for these solutions. In the left-hand figure we show the predicted regions of the masses of 
the lightest neutralino the lightest neutral Higgs scalar h° and the new neutral gauge 
boson Z 2 . In the right-hand figure we show the predicted cross section for the Z 2 decay 
into the dilepton final state (e + e~, // + //~) at the LHC with y/s = 14 TeV. The solutions 
are found only for more restricted values of 9 than those shown in Fig. 1. Figure 2 suggests 
that different annihilation processes of the lightest neutralino play an important role to 
realize its appropriate relic density for the CDM in different regions of 9. The variety of 
the predicted neutralino mass shows this aspect. 

As interesting results for mz 2 , we find that there exist upper bounds for mz 2 - Since 
the models considered in the present analysis have small tan f3 values, the mass bounds of 
the lightest neutral Higgs scalar can give a constraint on A as found from eq. (13). This 
constraint restricts allowed regions of u through the relation /i = Am. On the other hand, 
we can find the solutions only for restricted regions of /j, around 650 - 800 GeV in the 
present study. The upper bounds of mz 2 appear through these backgrounds. We also 
find that the lower bound of m Z2 can be less than 600 GeV in the allowed regions of 9. 
Since the singlino domination of the lightest neutralino makes the decay width of the Z 2 
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Fig. 2 The left-hand figure shows the masses of the lightest neutral Higgs, the lightest neutralino and the 
Zi boson in a GeV unit for various values of 9. The right-hand figure shows an expected cross section in 
a fb unit for the Z-i decay into the dilepton final states at the LHC with y/s = 14 TeV. These are derived 
for the parameters which satisfy all conditions discussed in the text. 



into the neutralino sector larger, the decay width into the dilepton final state becomes 
relatively smaller. Thus the constraint imposed by the analysis of the Z 2 search at the 
CDF seems to be escapable even for a smaller mz 2 than the usually discussed values. 
The predicted value of aB shows that Z 2 in the present models is easily detectable at the 
LHC. Combined searches of the Z 2 and the lightest neutral Higgs scalar will be useful to 
discriminate the models from other candidates beyond the SM. 



3.3 Details of annihilation processes of the lightest neutralino 

As mentioned in the previous part, different annihilation processes seem to play a crucial 
role for different regions of 6. To clarify these aspects, we examine the models defined 
by typical values of 9 in detail. We choose 6 = —0.4 and —0.253 as such examples. The 
features found in these examples are also seen in the models defined by other values of 8. 

In Fig. 3 we plot points in the (m^ 2 ,M^) plane, which satisfy all of the required 
conditions in the cases with /1 = 700 GeV and representative values of M^. The used 
values of are shown in the figures. In each figure we also show the regions surrounded 
by a dotted line (0 = —0.4) and a dashed line (0 = —0.253) where the lightest neutralino 
is dominated by the singlino S and also the conditions except for (29) are satisfied. The 
boundary expressed by a vertical line at a smaller m^ 2 is caused by the condition (28) 
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Fig. 3 Regions in the (mz 2 ,M s ) plane which satisfy all conditions discussed in the text, fi = 700 GeV 
is assumed. Values used for M^, and 9 are shown in each figure. 
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Fig. 4 Composition of the lightest neutralino in the case of 6 = —0.253, = 4.5 TeV and — 
200 GeV. 

which can give the lower bounds of mz 2 - O n the other hand, the boundary expressed by 
a vertical line at a larger mz 2 is caused by the mass bounds of the lightest neutral Higgs 
scalar given in (27). Since we fix a value of //, larger values of u make values of A smaller. 
As we can see from (13) and also as discussed in the previous part, smaller values of A can 
decrease the upper bounds of the Higgs mass and make it below the current experimental 
bounds at small values of tan f3. These features confine the mass of the new gauge boson 
into narrow regions as shown in each figure. Solutions for other values of fi have the 
similar features to those shown in this figure except that the solutions for a larger fi value 
move to the region with larger values of mz 2 - 

From Figs. 3 we find that all solutions are found in the regions where the lightest 
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Fig. 5 Masses of the lightest neutral Higgs and the lightest neutralino in the cases of 6 = —0.4 and 
—0.253. Myj is shown in each figure. 

neutralino is dominated by the singlino. In order to see the change of the composition 
of the lightest neutralino at the lower boundary of the regions of singlino domination, 
in Fig. 4 we plot values of \U^\ 2 defined in eq. (9) for mz 2 in the case of 9 — —0.253, 
M% = 4.5 TeV, and M-^ = 200 GeV. It shows that the singlino domination of the lightest 
neutralino suddenly turns into the bino domination within the narrow region of m^ 2 . The 
other cases with different values of and fi also show the similar behavior to this case 
on the same boundary. For larger values of M^y than 400 GeV we find the allowed regions 
in the same places in the {mz 2 , Mr) plane since the composition of the lightest neutralino 
is almost fixed there. 

In Fig. 5 we plot the masses of the lightest neutralino and the lightest neutral Higgs 
scalar obtained from the solutions shown in Fig. 3. Since larger values of Mj produce 
smaller masses for the singlino dominated lightest neutralino as discussed in section 2, 
we can find the correspondence of the solutions in Fig. 3 to the mass m^o of the the 
lightest neutralino shown in Fig. 5. We can find from Fig. 5 that there are three types 
of possibilities for the annihilation of the lightest neutralinos to explain the CDM abun- 
dance. They are characterized by the lightest neutralino mass. The first possibility is 
characterized by the relatively small mass such as ~ 45 GeV. It appears commonly in the 
all figures in case of 9 = —0.4. However, we cannot find this type of solutions in case of 
9 = —0.253. In these solutions the effective annihilation of the lightest neutralinos seems 
to be induced through a Z\ exchange. It can realize the appropriate CDM abundance 
only if m^o ~ rriz 1 /2 and the Higgsino components are suitably contained in the singlino 
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Fig. 6 Dominant processes which realize the suitable relic abundance fl x h 2 in the case of 6 = —0.4 and 
and Myy — 400 GeV. In order to realize three types of solution Ms is chosen such as 10.3 TeV for the Z\ 
exchange, 7.7 TeV for the Higgs exchange, and 1.9 TeV for the Z 2 exchange. 

dominated lightest neutralino. This conditions seems to be satisfied only in the restricted 
values of 9 such as —0.46 < 9 < — 0.34 as found in the left-hand figure of Fig. 2. The sec- 
ond one is characterized by the mass tuned into the narrow region around m/jo/2. These 
solutions appear as a result of the Higgs pole enhancement of the annihilation. The third 
one is characterized by the large neutralino masses such as ~ 170 — 180 GeV. In this case 
the annihilation is considered to be mainly induced by a Z 2 exchange. The large mass 
of the singlino dominated neutralino can enhance the factor bf in eq. (20) to realize the 
appropriate abundance. This type of solutions appear only for large values of such 
as 400 GeV or more. This is because the large is necessary for the heavy singlino 
dominated neutralino to be the lightest one. 

In order to check these points, we adopt the solutions for 9 = —0.4 and Myy = 400 GeV 
in Fig. 3 and calculate Q x h 2 by taking account of only the annihilation process discussed 
above for each case. We choose the values of M% as 10.3 TeV, 7.7 TeV and 1.9 TeV, where 
the three types of solutions appear. We plot the results in Fig. 6. This figure justifies the 
above discussion. Dominant components of the lightest neutralino in each case are found 
to be 



Z x exchange (Afe = 10.3 TeV) : W 6 -97.6%, jV 5 ~ 1.4%, 
h° exchange (M £ = 7.7 TeV) : W 6 -97.4%, Af 5 ~ 1.4%, 
Z 2 exchange {M$ = 1.9 TeV) : W 6 -91%, M x -7%, 
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Fig. 7 An expected cross section for the Zi decay into the dilepton final states at the LHC with 
y/s = 14 TeV in the cases of 6 = —0.4 and —0.253. /i = 700 GeV is assumed and My, is shown in each 
figure. 



For the heavy neutralino solutions the lightest neutralino can be heavier than the ordinary 
Z\. Then we may need to take account of additional new final states such as W + W~ and 
Z\h° which can be mediated by the Z 2 exchange. However, since these processes are gen- 
erally suppressed in comparison with the ones with final states // in the present models, 
these effects seem to be safely neglected in the annihilation of the singlino dominated 
lightest neutralino. Thus, the present results are considered to be a good approximation 
even in the case of m^o > mz 1 ■ 

In Fig. 7 we plot the predicted cross section aB of the Z 2 decay into the dilepton final 
states (e + e~, at the LHC with yfs = 14 TeV for the models defined by 9 = —0.4 
and —0.253. This figure also shows that the LHC can easily find Z2 in the present models. 
The models defined by the different values of 9 can also be distinguished from each other. 
In the previous part we discussed that the dilepton channel may be suppressed compared 
with the case where the lightest neutralino is composed of only the MSSM contents, 
since the Z2 decay into the neutralino pairs can be enhanced in the case of the singlino 
dominated lightest neutralino. This may be generally expected if S has the larger charge 
of U(l)a; compared with the charged leptons. In the present models this is the case and the 
aspect is shown in the figure for = 400 GeV in Fig. 7. In the case of = 400 GeV 
the appearance of new branches can be understood from this reasoning. They correspond 
to the solutions with the heavier lightest neutralino. Thus the Z 2 decay to the neutralino 
sector is suppressed and the cross section to the dilepton final states is enhanced. 
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4 Summary 



We have studied the phenomenology of the models with an extra U(l) which can give an 
elegant weak scale solution for the \i problem and have the singlino dominated lightest 
neutralino. Such models can be constructed by introducing an SM singlet chiral superfield 
S to the MSSM. The vacuum expectation value u of the scalar component of S generates 
the \x term and also breaks the extra U(l) symmetry. The singlino dominated lightest 
neutralino is realized by assuming the extra U(l) gaugino to be heavier than other gaug- 
inos. Under this assumption, the lightest neutralino can be dominated by the fermionic 
component of S even if the vacuum expectation value u is large enough to make the extra 
U(l) gauge field Z 2 satisfy the current experimental constraints. We studied the case 
with tan [3 = \jQi/Q2 where the ZZ' mixing constraint is automatically satisfied even for 
rather small value of u. Although this fact is known, the study of that case seems not to 
be done in the Z' phenomenology. 

The neutralino sector of the models is different from the currently known similar type 
models. In the NMSSM and the nMSSM the singlino dominated lightest neutralino usually 
appears for the small values of u. In the models with an extra U(l) in which the masses of 
the gauginos are assumed to be the same order, the lightest neutralino is expected to be 
very similar to that of the MSSM because of a large value of u. On the other hand, in the 
present models the singlino dominated lightest neutralino can appear in a very different 
situation from these and this can make their phenomenology distinguishable from others. 

We have focussed our study on the extra U(l) models which are derived from E 6 . And 
we have discussed typical phenomenological consequences in the cases where the singlino 
dominated lightest neutralino has the suitable relic density as a CDM candidate. We have 
shown that there are non-negligible parameter regions which satisfy the currently known 
phenomenological constraints as long as the model satisfies tan (3 ~ \jQi/Q2- The CDM 
candidate has been shown to have different natures from that in both the MSSM and 
other type of singlet extensions of the MSSM. They may be examined in future collider 
experiments. We have also predicted the masses of the lightest neutral Higgs scalar, the 
new neutral gauge boson and the lightest neutralino. The predicted cross section for the 
Z 2 decay into the dilepton pairs shows that the Z 2 boson in the models can be easily 
found at the LHC. Searches of the neutral Higgs scalar and the Z 2 boson at the LHC are 
expected to give us useful informations for the models with an extra U(l) which may give 
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the solution for the \i problem. 
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Appendix 

In this appendix we present the interaction Lagrangian relevant to possible decay 
processes of Z 2 and also formulas for the decay width for them [14]. 

Interaction terms relevant to this decay process are 

<& = \ (Q(h) + Q(fn)) , ql = \ (Q(h) - QUr)) ■ (30) 



The decay width can be calculated as 



_ g 2 x m Z2 



JN 



mi , x 

4-f , 31 



where c/ equals to 1 for leptons and 3 for quarks. Because of various reasons, the extra 
fermions in the present models can be considered to be sufficiently heavy and then ne- 
glected in the final states for the Z 2 decay. In that case it is justified to confine / into 
ordinary quarks and leptons. 
(ii) Z 2 - f\ R f\ R 

Mass matrices of the sfermions have a similar structure to that in the MSSM except that 
there are additional contributions from the extra U(l) D-term. It can play a crucial role 
since it may give large negative contributions depending on the models. They can be 
written as 



M) 



JLJL JLJR 

, Mf f Ml -. 

\ JLJR JrJr 



(32) 



where each component of this matrix is a 3 x 3 matrix with respect to flavors. However, 
if flavor mixing is small enough, each flavor can be treated independently and M. 2 j can be 
reduced into a 3 x 3 matrix. In that case the component of .Mj can be expressed as 

M hf L =fn l + l m /! 2 + ( T 3 - QfSw) m z cos 2 P + Qx{fh)rri 2 D , 

M )r]r =^o + I™/| 2 + Q f s 2 v m 2 z cos2P + Q x (f R )m 2 D , 

M} Lf ~ R = -m* f (A*-\uR f ), (33) 

where (T 3 ,Rf) takes (|,cot/3) and (— |,tan/3) for an up- and down-sector of squarks and 
sleptons, respectively, to/ and Qf are the mass and the electric charge of a fermion /. 
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We assume the universality of soft scalar masses mj and A parameters. The D-term 
contribution of the extra U(l) has the expression 

™r> = \ffc {Qivl + Qiv\ + Qsu 2 ) . (34) 

We should note that this contribution may be negative to realize the light sfermions 
depending on the value of introduced in eq. (24) which defines the extra U(l) models. 
Interaction terms relevant to this decay process are given by 

C = 9 -^(ql ± ql)~fl R i d, f LtR Z$, (35) 



where Vf and a/ are defined in eq.(30). The decay width of this process can be written as 
r- -, = c ^^ q f ±q fy(i-4 r ^-) . (36) 

In this derivation the LR mixing Mj j in eq. (32) is neglected. This is considered to be 
a good approximation except for the stop sector. 

(iii) Z 2 -> H+H- 

The charged Higgs sector has the same structure as that of the MSSM except that there 
are additional contributions to the mass eigenvalue generated by the coupling \SHiH 2 as 
in the NMSSM. The mass eigenvalue and its eigenstate can be expressed as 

ml ± =m ^(i-^p\+ tf± = H f sin (3 + Hf cos/3. (37) 

V 9i J sin 2/3 

The interaction Lagrangian relevant to this process is given by 

C = ^ (Q l sin 2 (3 - Q 2 cos 2 (3) H+i H-Z%, (38) 

where G ± are would-be Goldstone bosons. Using this interaction Lagrangian, we can 
derive the decay width for this process as 

2 / 2 \ 3 / 2 

T H+H - = 9 -f^(Qi sin 2 (3-Q 2 cos 2 1 - 4^ . (39) 

(iv) Z 2 -> W ± H^ 

The relevant interaction Lagrangian for this process is 



£ = y (Qi + Q2) sin/9 cos (3 



(m w H-W+ + H~i <9 M G + )Z% + h.c. 
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(40) 



Using this Lagrangian, the decay width for this process can be calculated as 
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(Qi + Q 2 ) 2 sin 2 /3cos 2 /? 



^2 



1/2 



X 



1 + 2 5m 2 w -m 2 H± + (m^ - m 2 H± f 



m 



z 2 



m 



z 2 



(41) 



(v) Z 2 -> 

The relevant interaction Lagrangian is written as 

£ = £ 2^m z (Q 1 lQ cos/? - Q 2 2a sin/3)^ + d^G°)(j) a Z%. 

a=l 1 



(42) 



where we denote a would-be Goldstone boson as G°. The matrix (9 and the mass eigen- 
values are defined in eq. (12) and also in the text below it. The decay width for this 
process can be derived as 



Z(j> a 



3 2 

E 9 -j~z{QiO Xa cos /? - QA sin pf 



1 - 2 



m 



z 2 



rn 



z 2 



1/2 



X 



(43) 



a=l 

m 2 Z2 m% 2 

(vi) Z 2 ->• a P4 

Pa is the CP odd Higgs scalar. The relevant interaction Lagrangian for this process is 
written as 

C=±%[ QiOj^- + ft©*^ + Q S Os«^^) 2f i>x 5 A., (44) 



o=l 



N ^« -* u ^ N 

where N is defined in eq. (16) and in the text below it. Using this, we can derive the 
decay width as 

g 2 x m Z2 1 



4 48tt N 2 



(Qi0 la u sin (3 + Q 2 2a u cos (3 + Q s 3a v sin (3 cos f3Y 



x 



m 2 , + mi 
1-2 " Pa + 
ml 



*Z 2 m Z 2 

(vii) Z 2 -> 

The neutralino sector is discussed in detail in the text. The interaction Lagrangian for 
the neutralinos x° a an d Z 2 is given by 



3/2 



(45) 



a,b=l 



(46) 
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where an effective coupling g ab is defined as 



got = E 

j=4 Z 



(47) 



The definition of the matrix [/ is given in eq. (9). Using this interaction Lagrangian, the 
decay width is derived as 



9ab m Z 2 

12tt 



1 - 2*o* 



m a + m 2 {m a — m 2 ) 3m a m b 



Z 2 



X 



i _ 2 mg + mg + (m 2 - m 2 b ) 2 



«z 2 ™z 2 

where m a stands for the mass eigenvalue of the neutralino x° a - 
(viii) Z 2 -> C+C^ 

The chargino sector has the same structure as that of the MSSM. The mass terms can be 
expressed as 



1/2 



(48) 



-i\ w -, H 1 



Myy \f2m-w sin (3 \ I —iX 



w+ 



\[2m w cos (3 



-Xu 



+ h.c. 



(49) 



If we define the mass eigenstates as 




V + 



V- 



-i\\v+ 
Hi 

-iX w - 
Hi 



U_ = 



COS 


— sin 4> + 


sin 


COS (j) + 


COS(/>_ 


— sin 0_ 


sin 0_ 


COS(/>_ 



(50) 



the mixing angles (f>± have analytic expressions 

2^/2m w (Myy sin (3 — Xu cos (3) 



tan20 4 
tan20_ 



2 — M 2 ' 



2771^ (sin 2 /3 - cos 2 (3) + A 2 m 
2y / 2mvy(M T y cos (3 — Xu sin /3) 



2m 2 K (cos 2 /3 - sin 2 (3) + X 2 u 2 - M\ 
Using this basis, the mass terms can be transformed into 



(51) 




+ h.c, 



(52) 
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where the mass eigenvalues are represented as 



mi = (Myy sin 0_ + V2mw cos (3 cos 0_) cos + 

— ( \^2mw sin /9 sin 0_ — \u cos 0_ ) sin + , 
m 2 = (M^ cos 0_ — v^mvK cos (5 sin 0_) sin + 

+ ( V2m w sin /3 cos 0_ + Am sin 0_ ) cos + . 



(53) 



Since only the Higgsinos H^ 2 interact with Z 2) the interaction Lagrangian relevant to the 
Z 2 decay is given by 



£ = TT 12 C aln{Vab + aablb)C b Z^ 



(54) 



a,b=l 



where = (ip^ ,x„), and effective couplings v ab and a ab are defined as 
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■ 2 i . Q2 . 2 i 

sin 0_ + — sin 0_| 
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V22 
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2 1 , Q2 2 1 
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Vl2 
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= 0,21 = 


sin 0_ cos 0_ 
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^ ■ A 

sin 
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Qi . 2 

— sin < 
2 

Qi 2 

■ — cos 
2 



Q2 . 2 
— sin 1 
2 

Q2 5 

cos 

2 



(55) 



Using these couplings, the decay width for this process is derived as 



C a Cb 



9x m z 2 
4 48tt 



ml + m 2 b (ml-ml) 2 ^ 



2 m a m b 



y ab + -ib) ( 1 - - v v 4 b; ) + 3(«i - < br -,. 



m 



z 2 J 



X 



j _ n m l + ^ + (™a ~ ^b) 2 



ffl; 



^ 2 m z 2 

Finally, the branching ratio B(Z 2 — > X7) used in the text is defined by 

r(z 2 -> xf) 



1/2 



(56) 



5(Z 2 -> X7) 
where the total decay width r tot of Z 2 is 



(57) 



tot 



r tot - 12 T ff + 12 T fL,Rft R + v h+h- + r w ± H ^ 
/ / 

3 3 6 2 

+ 2 ^ + 12 T <PjPA + 12 T NaN b + 12 V C a C b - 

i=l j=l a,b=l a,b=l 



(58) 
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